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The coronavirus family of positive-strand RNA viruses includes important pathogens of livestock, compan- 
ion animals, and humans, including the severe acute respiratory syndrome coronavirus that was responsible 
for a worldwide outbreak in 2003. The unusually complex coronavirus replicase/transcriptase is comprised of 
15 or 16 virus-specific subunits that are autoproteolytically derived from two large polyproteins. In line with 
bioinformatics predictions, we now show that feline coronavirus (FCoV) nonstructural protein 16 (nsp16) 
possesses an S-adenosyl-L-methionine (AdoMet)-dependent RNA (nucleoside-2’O)-methyltransferase (2’O- 
MTase) activity that is capable of cap-1 formation. Purified recombinant FCoV nsp16 selectively binds to short 
capped RNAs. Remarkably, an N7-methyl guanosine cap ("“°GpppAC,.,) is a prerequisite for binding. 
High-performance liquid chromatography analysis demonstrated that nsp16 mediates methyl transfer from 
AdoMet to the 2’0 position of the first transcribed nucleotide, thus converting 7“°GpppAC;., into 
7™eGpppAsomeC3.¢: The characterization of 11 nsp16 mutants supported the previous identification of resi- 
dues K45, D129, K169, and E202 as the putative K-D-K-E catalytic tetrad of the enzyme. Furthermore, residues 
Y29 and F173 of FCoV nsp16, which may be the functional counterparts of aromatic residues involved in 
substrate recognition by the vaccinia virus MTase VP39, were found to be essential for both substrate binding 
and 2'0-MTase activity. Finally, the weak inhibition profile of different AdoMet analogues indicates that nsp16 
has evolved an atypical AdoMet binding site. Our results suggest that coronavirus mRNA carries a cap-1, onto 
which 2’0 methylation follows an order of events in which 2’O-methyl transfer must be preceded by guanine 


N7 methylation, with the latter step being performed by a yet-unknown N7-specific MTase. 


The family Coronaviridae, comprising the genera Coronavi- 
rus and Torovirus, belongs to the order Nidovirales, a lineage of 
positive-strand RNA viruses that also includes the Arteriviridae 
and Roniviridae families (for a review, see reference 28). Coro- 
naviruses (CoV) are frequently associated with respiratory and 
enteric diseases in humans, livestock, and companion animals. 
In recent years, they received worldwide attention following 
both the 2003 outbreak caused by the emerging severe acute 
respiratory syndrome CoV (SARS-CoV) (17, 43, 57) and the 
subsequent identification of several other novel family mem- 
bers, including two additional human pathogens (63). On the 
basis of antigenic and genetic analyses, CoV have been divided 
into three groups (29), with group 3 representing avian CoV 
and group 2 including viruses like murine hepatitis virus 
(MHV), the human CoV (HCoV) OC43 and HKU1, and 
SARS-CoV. CoV group 1 includes, among others, the HCoV 
229E and NL63, the porcine transmissible gastroenteritis virus, 
and canine and feline CoV (FCoV). Infection with FCoV is 
very common in cats and kittens. Although FCoV infection 
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usually is transient, resulting in a mild gastrointestinal disease, 
the virus may persist in a significant percentage of cases (1). A 
small proportion of FCoV-infected cats develop a lethal, im- 
mune-mediated disease that is known as feline infectious peri- 
tonitis (FIP). The basis for virulence has remained controver- 
sial, with one suggestion being that FIP virus (FIPV) arises 
when a persistently infecting, enteric FCoV acquires mutations 
that increase its virulence (81). 

As for other nidoviruses, the large genome of CoV (27 to 32 
kb) is polycistronic, with about two-thirds being occupied by 
two large replicase open reading frames (ORFla and ORF 1b) 
that encode the viral nonstructural proteins (nsps). The genes 
downstream of the replicase ORFs encode structural and vi- 
rus-specific accessory proteins (for a review, see reference 28). 
Genome expression starts with the translation of ORFla and 
ORF 1b, presumably by a cap-dependent mechanism (44, 80), 
with the expression of ORF 1b involving a —1 ribosomal frame- 
shift (9). The two resulting replicase polyproteins, ppla and 
pplab, are processed by two or three viral proteases to gener- 
ate 16 end products, termed nsp1 to nsp16 (28, 76, 87). These 
cleavage products assemble into a large, membrane-anchored 
multienzyme complex, termed the replication-transcription 
complex, that mediates all functions required for genome rep- 
lication and subgenomic mRNA synthesis (38, 78, 79). The 
replication-transcription complex includes enzyme functions 
commonly found in positive-strand RNA viruses, like an RNA- 
dependent RNA polymerase (nsp12 [13]), proteases (nsp3 and 
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nsp5 [4, 49]), and a helicase/RNA triphosphatase (nsp13 [40, 
73]). Additionally, the CoV genome encodes a set of RNA- 
processing activities that either are unique to certain nidovirus 
subgroups or are found in only a few other groups of RNA 
viruses (76). These functions include an ADP-ribose phosphatase 
(X or macro domain in nsp3 [62]), a recently discovered putative 
RNA primase (nsp8 [37]), an exoribonuclease (nsp14 [54]), and a 
nidovirus uridylate-specific endoribonuclease (NendoU in nsp15 
[39)). 

This group of enzymes also includes the functionally unchar- 
acterized nsp16, which was predicted to be an S-adenosyl-L- 
methionine (AdoMet)-dependent RNA (nucleoside-2'O)- 
methyltransferase (2’O-MTase) (53, 76, 82). It contains a 
highly conserved catalytic tetrad (K-D-K-E) that is a hallmark 
of RNA 2'O-MTases (12, 19, 22, 23). A three-dimensional 
model of the MTase core of SARS-CoV nsp16 was generated 
by a structure prediction server (3D jury meta predictor) using 
the 2'O-MTase domain of the reovirus protein \2 as a tem- 
plate (82). Based on the role of its homologues in other RNA 
viruses, nsp16 was postulated to be involved in mRNA capping 
(76, 82). Alternatively, but not mutually exclusively, it was 
proposed that nsp16 2’O methylates selected nucleosides in 
virus and/or cellular RNAs as part of a pathway that could 
involve other unique CoV RNA-processing enzymes (76). 
Somewhat compatible with this hypothesis, RNA cleavage by 
the nsp15 NendoU was found to be inhibited by 2’O-methyl- 
ation (39). Whereas the exact role of nsp16 during CoV rep- 
lication still is unknown, its functional importance was sup- 
ported by mutagenesis experiments using a SARS-CoV 
replicon system (3). The deletion of the nsp16 coding sequence 
blocked RNA synthesis, whereas a single mutation in the cat- 
alytic tetrad reduced replicon-driven mRNA synthesis to about 
10% of the level for the wild type (wt). Moreover, the pheno- 
type of temperature-sensitive MHV mutants suggested that 
nsp16 has an essential role in either the synthesis or stability of 
viral RNA or in controlling a cellular function that is able to 
limit virus replication (70). Elucidating the specific biochemi- 
cal properties of nsp16 clearly would help to assess its role and 
importance in the viral life cycle. 

The RNA cap is a unique structure found at the 5’ end of 
eukaryotic cellular and many viral messenger RNAs. In eu- 
karyotic cells, the cap structure protects mRNA from degra- 
dation by 5’ exoribonucleases and enhances the initiation of 
mRNA translation (24, 75). For nascent cellular transcripts, 
the addition of the cap-0 structure is a cotranscriptional event 
that occurs in the nucleus. Cap-0 formation generally requires 
three sequential enzymatic activities: (i) an RNA 5’-triphos- 
phatase (RTPase) that removes the 5’ y-phosphate group of 
the mRNA; (ii) a guanylyltransferase (GTase), or capping en- 
zyme, that catalyzes the transfer of GMP to the remaining 
5'-diphosphate end; and (iii) an AdoMet-dependent (guanine- 
N7)-MTase (N7-MTase) that methylates the cap at the N7 
position. Whereas lower eukaryotes, including yeast, employ a 
cap-0 structure, higher eukaryotes convert cap-0 into cap-1 or 
cap-2 structures (46) by means of a nuclear AdoMet-depen- 
dent 2'O-MTase that methylates the ribose 2’O position of the 
first and second nucleotide of the mRNA, respectively. 

Many viruses that replicate in the cytoplasm encode their 
own RNA capping machinery. Some of these viruses, such as 
the positive-sense single-stranded RNA (ssRNA) flaviviruses 
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(20, 66) and DNA poxviruses (84), seem to have adopted the 
sequential four-step mechanism used in eukaryotic mRNA 
cap-1 formation. However, the molecular and genetic organi- 
zation of the enzymatic activities involved in RNA capping 
varies between virus groups. For example, in the poxvirus vac- 
cinia virus (VV), cap-0 formation is catalyzed by a single 95- 
kDa protein encoded by viral gene D1 (21, 31, 74, 84). The 
subsequent methylation of capped RNA at the 2’O position 
requires VP39, a bifunctional protein that also directs 3’ poly- 
adenylation (71). In contrast, in flaviviruses, an RTPase activity 
is found in the C-terminal domain of the multifunctional he- 
licase protein NS3 (7, 8), whereas the two MTase activities (the 
N7- and 2'O-MTases) reside in the N-terminal domain of the 
RNA-dependent RNA polymerase subunit NS5 (19, 66). 
There also are RNA viruses with capping mechanisms that 
deviate dramatically from the canonical pathway. Alphaviruses 
may use an alternative pathway for mRNA capping, in which 
the GTP molecule is methylated before being transferred to 
the 5’-diphosphate end of viral RNAs (2). Yet another uncon- 
ventional mechanism is employed by the rhabdovirus vesicular 
stomatitis virus (VSV), which involves a unique polyribonucleo- 
tidyltransferase activity to transfer the monophosphate mRNA 
onto GDP derived from GTP (55). 

CoV not only have a genomic mRNA but also produce an 
extensive nested set of subgenomic mRNAs, a property that 
emphasizes the importance of the presumed RNA capping 
process (for recent reviews, see references 56 and 69). On the 
basis of T1-oligonucleotide fingerprinting, Lai and Stohlman 
previously claimed that the 5’ end of MHV mRNAs carries a 
cap structure (44). Moreover, using a cap-specific monoclonal 
antibody and exoribonuclease protection assays, it was dem- 
onstrated that both genomic and subgenomic mRNAs of 
equine torovirus carry a 5’ cap (80). However, the cap struc- 
ture and the enzymes involved in CoV RNA capping and their 
mechanisms of action have not yet been characterized. The 
multifunctional nsp13 helicase subunit possibly is involved, 
which was previously shown to carry an RNA 5’-triphosphatase 
activity in the case of both group 1 and group 2 CoV (HCoV- 
229E and SARS-CoV) (40, 41). As for several other groups of 
RNA viruses, the GTase involved has not yet been identified. 
In addition to the predicted nsp16-mediated 2’O-MTase activ- 
ity (76, 82), it was suggested that the SUD (for SARS-CoV 
unique domain) of SARS-CoV nsp3 exhibits N7-MTase activ- 
ity (26). This prediction, which has not been verified experi- 
mentally, is weakened by the fact that this domain is conserved 
only in group 2b CoV (A. E. Gorbalenya, unpublished data). 
Thus, the identification of the N7-MTase as well as other CoV 
capping factors remains to be achieved. 

In this paper, we provide the first experimental evidence for 
the AdoMet-dependent 2'O-MTase activity of a CoV nsp16. 
Using recombinant FCoV nsp16 and short capped and un- 
capped RNAs, we show that this enzyme has selective RNA 
binding properties and demonstrate that it is a cap-0 ("“°Gp 
ppAC;) binding protein. Using site-directed mutagenesis, we 
confirm the essential role of the conserved K-D-K-E catalytic 
tetrad for mRNA cap 2’O-MTase activity and identify aro- 
matic residues in the nsp16 N terminus that may play a key role 
in cap-0 binding. Finally, the weak inhibition profile of 
AdoMet analogues that are known potent MTase inhibitors 
suggests the presence of a unique methyl donor binding site in 
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an otherwise conserved three-dimensional MTase domain or- 
ganization. 


MATERIALS AND METHODS 


Cloning of the FCoV nsp16 gene. FCoV strain FIPV WSU-79/1146 (GenBank 
accession no. DQ010921) (18) was kindly provided by Stuart Siddell (University 
of Bristol, United Kingdom) and was used to infect Crandell-Reese feline kidney 
cells. Intracellular RNA was isolated from infected cells and was used for the 
reverse transcription-PCR amplification of the nsp16 coding sequence (genome 
residues 19307 to 20206). The PCR primers used (sense, 5’°-GGGGACAAGTT 
TGTACAAAAAAGCAGGCTTCGAAGGAGATAGAACCATGAAACATC 
ACCATCACCATCACAGTTTAGAAAATGTGGCTTATA-3’; and antisense, 
5'-GGGGACCACTTTGTACAAGAAAGCTGGGTCTTATTGTAGTTTTG 
GGTAGAAGGTT-3’) specified an N-terminal, in-frame hexahistidine tag and 
recombination sequences for use in the Gateway cloning system (Invitrogen). 
Following cloning into entry vector pDONR201 and sequence verification, the 
nsp16 gene was transferred to Escherichia coli expression vector pDest14 (In- 
vitrogen) to give pDest14/6His-nsp16. Nsp16 mutagenesis was performed using 
the QuickChange site-directed mutagenesis kit (Stratagene) according to the 
manufacturer’s instructions, followed by sequence verification (Millegen, 
France). 

Reagents. AdoMet and adenosine-homocysteine (AdoHcy) were purchased 
from New England Biolabs. AdoMet analogues were obtained from Sigma- 
Aldrich. They were dissolved in dimethylsulfoxide (DMSO) and stored as 10 mM 
stock solutions at —20°C. «-3*P-labeled cytosine 5’-triphosphate (3,000 Ci/ 
mmol), uniformly labeled H]GTP (5.20 Ci/mmol), and 5’ triphosphate were 
purchased from GE Healthcare. 

Expression and purification of NSSMTasepy, hMTase, and FCoV nsp16. 
Recombinant NSSMTasepy, corresponding to residues 1 to 296 of dengue virus 
(DV) NSS, was expressed and purified as described previously (19). The Human 
N7-guanine MTase (hN7-MTase) cDNA was a kind gift from Aaron J. Shatkin 
(Center of Advanced Biotechnology and Medicine, Piscataway, NJ). It was pro- 
duced and purified as described previously (58). 

E. coli C41(DE3) (Avidis SA, France), transformed with the pLysS plasmid 
(Novagen), was transformed with pDest14/6His-nsp16 and grown in Terrific 
broth (TB) containing ampicillin and chloramphenicol. When the optical density 
at 600 nm reached 0.6, isopropyl-B-b-thiogalactopyranoside was added to a final 
concentration of 0.5 mM, and expression was allowed to proceed for 18 h at 17°C. 
The bacterial cell pellet was resuspended in lysis buffer (50 mM Tris-HCl, pH 8.0, 
300 mM NaCl, 5% glycerol, antiprotease cocktail [complete; Roche]) supple- 
mented with 10 mM imidazole, 100 jg/ml lysozyme, 0.25 g/ml Dnase I, and 
0.1% Triton X-100. After lysis by sonication and clarification, immobilized metal 
affinity chromatography was used for the first purification step (chelating Sepha- 
rose fast-flow resin; GE Healthcare). FCoV nsp16 was eluted with 250 mM of 
imidazole in lysis buffer. Fractions containing nsp16 then were filtered (0.45-~m 
filter), loaded onto a Hi Load 16/60 Superdex 200 gel filtration column (GE 
Healthcare), and eluted with 10 mM HEPES, pH 7.5, 300 mM NaCl, 2 mM 
dithiothreitol (DTT), 50 mM arginine, and 50 mM glutamate. The fractions 
containing nsp16 were concentrated to 2 mg/ml on Amicon ultra 15 centrifugal 
filter units (10 kDa; Millipore) and stored at —20°C in the same buffer supple- 
mented with 50% glycerol. The identity of wt nsp16 and mutants was checked by 
matrix-assisted laser desorption ionization—time of flight mass spectrometry after 
trypsin digestion (with wt nsp16 showing 43.2% of sequence coverage), by dy- 
namic light scattering, and by sodium dodecyl sulfate-polyacrylamide gel elec- 
trophoresis (SDS-PAGE). All mutants, except for the poorly expressed W175A, 
were purified by the two-step procedure described above for wt nsp16. Mutant 
proteins showed gel filtration elution profiles similar to those of the wt control 
and were recovered as dimers from the column. Comparable amounts of each 
mutant protein were analyzed by SDS-PAGE (see Fig. 7B). Their electro- 
phoretic mobility patterns were identical, except for that of the D129A mutant, 
which showed slightly aberrant migration in the gel (a 1- to 2-kDa decrease in 
apparent molecular mass). 

RNA synthesis and purification. Capped and noncapped RNAs (7°Gpp 
pAC;, GpppACs, and pppACs) were synthesized in vitro using bacteriophage T7 
DNA primase and were purified by high-performance liquid chromatography 
(HPLC) as described previously (58). **P-labeled pppAC,, ™“°GpppAC,,, or 
GpppAC,, RNA was synthesized as described previously (20). Reactions were 
stopped by adding RNase-free DNase I (500 U/ml for 30 min at 37°C; GE 
Healthcare) and proteinase K (100 g/ml for 16 h at 37°C; Invitrogen) to remove 
the DNA template and the primase present in the reaction mixture. Proteinase 
K was inactivated at 70°C for 5 min. After 15 min of centrifugation at 16,000 x 
g in order to remove insoluble material, RNAs were resuspended in 100 pl of 10 
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mM Tris-HCl, pH 7.5, 50 mM NaCl, 2.5% glycerol, 500 g/ml bovine serum 
albumin (BSA) and stored at —20°C. 

MTase activity assays. Reactions for MTase activity assays were performed in 
40 mM Tris-HCl, pH 7.5, 5 mM DTT, 2 4M ™*°GpppAC, or GpppAC,, 5 uM 
AdoMet, and 0.03 mCi/ml [*H]AdoMet (GE Healthcare). MgCl, was added at a 
concentration of 1 mM for the standard nsp16 MTase assay. The reducing agent 
DTT (1 mM) was found to stabilize nsp16 activity over time (data not shown). 
NSSMTasepy, hN7-MTase, and FCoV nsp16 were added at concentrations of 
500 nM, 200 nM, and 3 uM, respectively, and reaction mixtures were incubated 
at 30°C and stopped after 4 h by a 10-fold dilution of the reaction mixture in 100 
pM ice-cold AdoHcy. Samples were kept on ice and then transferred to glass- 
fiber filtermats (DEAE filtermat; Wallac) by a Filtermat Harvester (Packard 
Instruments). Filtermats were washed twice with 0.01 M ammonium formate, pH 
8.0, twice with water, and once with ethanol, dried, and transferred into sample 
bags. Betaplate Scint (Wallac) scintillation fluid was added, and the methylation 
of RNA substrates was measured in counts per minute (cpm) by using a Wallac 
1450 MicroBeta TriLux liquid scintillation counter. For inhibition assays, we set 
up the reaction as described above with a short RNA substrate ("“°GpppAC;) 
in the presence of various concentrations of candidate inhibitors. Enzymes and 
RNA substrates were mixed with the inhibitor before the addition of AdoMet to 
start the reaction. The final concentration of DMSO in the reaction mixtures was 
below 5%, and control reaction mixtures without inhibitor contained corre- 
sponding DMSO concentrations. Reaction mixtures were incubated at 30°C for 
4h. The samples then were transferred to glass-fiber filtermats (DEAE filtermat; 
Wallac) and counted as described above. The ICs, (inhibitor concentration at 
50% activity) value of GTP, 7“°GTP, and AdoHcy were determined using 
Kaleidagraph. Data were adjusted to a logistic dose-response function, % activity = 
100/(1 + [IJ/IC5o)’, where b corresponds to the slope factor that determines the 
slope of the curve and [I] corresponds to the inhibitor concentration (15). 

For HPLC analysis, reactions were performed in the absence of [*H]AdoMet, 
and samples were mixed with 200 pl of triethylammonium bicarbonate solution 
(TEAB) (0.05 M) before HPLC analysis on a Waters model 600 gradient HPLC 
system. The column assembly consisted of a precolumn (Delta-pak C,,; 100 A, 
5 ym, 3.9 < 20 mm) and a separation column (Nova-pak C,s; 4 wm, 3.9 x 150 
mm). For the on-line cleaning procedure, both columns were installed in parallel 
on a 7000 Rheodyne two-valve system (Interchim). Eluent A was a 0.05 M 
solution of TEAB (pH 7.4), and eluent B was a 1:1 (vol/vol) mixture of aceto- 
nitrile (HPLC grade; Carlo Erba SDS, France) and TEAB (final concentration, 
0.05 M; pH 7.4). Separations were run at a flow rate of 1 ml/min and started with 
a 5-min elution (100% eluent A) on the precolumn to remove proteic material. 
The gradient started after 5 min at 100% eluent A, with an increase to 10% 
eluent B after 25 min and to 30% eluent B after 45 min. 

Analysis of the 7”“°GpppAC; product by enzymatic digestion and HPLC. A 
preparative-scale reaction (800 jl) was performed using the standard conditions, 
and the reaction was stopped after overnight incubation. HPLC purification was 
performed as described above, and peaks corresponding to the reaction products 
were collected. Lyophilization was repeated three times, and 3 nmol of the 
reaction products was digested in a 60-pl reaction volume containing 50 mM 
Tris-HCl, pH 8.5, 5 mM MgCl, 0.3 U of nucleotide pyrophosphatase type II 
from Crotalus adamanteus (0.03 U in a 60-ul reaction volume; the enzyme 
preparation contained phosphodiesterase I side activity; Sigma) and 60 U of calf 
intestinal phosphatase (New England Biolabs). The digestion was allowed to 
proceed for 25 min at 37°C and was stopped by being heated to 95°C for 5 min. 
The digestion products then were separated by HPLC as described previously 
(58). 

MTase binding assays. Binding assays were performed for 2 h at 4°C in 
binding buffer (10 mM Tris-HCl, pH 7.5, 50 mM NaCl, 2.5% glycerol, and 500 
pg/ml BSA) in a total volume of 150 yl, as described previously (20). For 
competition experiments, incubation mixtures contained increasing concentra- 
tions of MgCl, 7“°GTP, GTP, GpppA, or ”“°GpppA. NSSMTasepy and FCoV 
nsp16 bound to nickel-nitrilotriacetic acid (Ni-NTA) beads were produced as 
described above. For each experiment, we used 30 yl of NSSMTasepy beads (~4 
ug/l) and 10 pl of 3?P-radiolabeled RNA. Accordingly, 1 nmol of RNA was 
incubated with beads containing approximately 4 nmol of protein. After three 
washes with binding buffer, the RNA bound to the beads was liberated using 10 
pl of formamide-EDTA gel loading buffer and separated by polyacrylamide gel 
electrophoresis (14% acrylamide-bisacrylamide [19:1], 7 M urea) in TTE buffer 
(89 mM Tris-HCl, pH 8.0, 28 mM taurine, 0.5 mM EDTA). RNA bands (with the 
input corresponding to 1/10 or 1/20 of the RNA incubated with NSSMTasepy, 
beads) were visualized using photostimulated plates (fluorescent image analyzer 
FLA3000; Fuji). 

Three-dimensional model of the nsp16 MTase core. Different Web servers 
employing threading and homology modeling (CPHmodels, Fugue, Phyre, 
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FIG. 1. Purification of FCoV nsp16 and AdoMet-dependent MTase assay. (A) A gel summarizing the purification of nsp16. After two 
chromatographic steps (see Materials and Methods), purified FCoV nsp16 protein was analyzed using 12% SDS-PAGE and the gel was stained 
using Coomassie blue. Lane 1 corresponds to the molecular mass marker and lane 2 to the final nsp16 protein preparation. (B) AdoMet-dependent 
MTase activity of hN7-MTase, NSSMTasepy, and FCoV nsp16. Equal amounts of the different enzymes were incubated with GpppAC; and 
™eGpppAC; in the presence of [*H]AdoMet. The methyl transfer to the RNA substrate was monitored during a 120-min time course and detected 
using a filter binding assay that measured the transferred amount of radioactivity in cpm. (C) Time course of the methylation reaction. FCoV nsp16 
was incubated with "“°GpppAC; in the presence of [*H]AdoMet. The methyl transfer to the RNA substrate was monitored during 360 min and 


detected as described above. 


ESyPred3D, HHpred, and LOOPP) were used to generate a model of FCoV 
nspl6. The closest template proposed by several servers was the rRNA 2'0O- 
MTase FtsJ (Protein Data Bank identity no. [PDB ID] 1EIZ [11]), which served 
as the basis to generate a model using the HHpred server (http://toolkit 
.tuebingen.mpg.de/hhpred). K45 was added manually using TURBO (68). 


RESULTS 


Purified recombinant FCoV nsp16 protein is an AdoMet- 
dependent MTase. The CoV nsp16 previously was predicted to 
be a 2'O-MTase, therefore the enzyme was postulated to par- 
ticipate in cap-1 ("“°GpppN2'O,,.) formation (76, 82). To 
experimentally verify the MTase activity, a cDNA sequence 
encoding FCoV nsp16 (300 amino acids; FCoV strain FIPV 
WSU-79/1146 [18]) was cloned and expressed in EF. coli to 
produce a recombinant nsp16 that was (His), tagged at its N 
terminus. The protein was purified by immobilized metal af- 
finity chromatography, followed by size-exclusion chromatog- 
raphy in the presence of arginine and glutamate (each at 50 
mM) in order to avoid protein precipitation (27). Upon gel 
filtration, the protein eluted as a single peak, corresponding to 
a size of about 65 kDa. In view of the calculated molecular 
mass of nsp16 and its migration as a single 35-kDa band during 


SDS-PAGE (Fig. 1A), this indicated that the purified protein 
probably was dimeric in solution. The identity of the recombi- 
nant protein was confirmed by matrix-assisted laser desorption 
ionization—time of flight mass spectrometry after trypsin diges- 
tion (data not shown). 

To evaluate whether FCoV nsp16 exhibits MTase activity, 
we first incubated the purified protein with different short 
capped or uncapped RNA oligonucleotides (’“°GpppACs, 
GpppAC.,, and pppAC;) in the presence of the radiolabeled 
methyl donor [*H]AdoMet. The first two nucleotides (AC) of 
the substrates were identical to the authentic 5’-terminal nu- 
cleotides of the FCoV genome (5'-ACUUUU. . .). As positive 
controls for MTase activity, we used hN7-MTase and 
NSSMTasepy (19, 58). Aliquots of the reaction mixtures were 
absorbed onto filters and washed, and the radioactivity that 
remained associated with the RNA substrate was measured to 
quantify its methylation. 

No label was transferred to the uncapped pppAC; substrate, 
using either nsp16 or both control MTases (data not shown). 
As expected (58), NSSMTase,y shows MTase activity on both 
capped RNA substrates, whereas hN7-MTase methylated Gp 
ppAC; but not ™“°GpppAC; (Fig. 1B). In contrast, FCoV 
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nsp16 did not catalyze methyl transfer to either the guanine N7 
or the 2’0 position of the GpppAC, substrate. However, it was 
active on the capped 7@*GpppAC, substrate and, thus, acted 
as an AdoMet-dependent mRNA cap MTase (Fig. 1B). From 
these observations, we conclude that FCoV nsp16 can transfer 
a methyl group to the 2’O position of the first and/or second 
nucleotide of an N7-methylated short RNA substrate. A time 
course experiment using FCoV nsp16 showed that, under our 
standard assay conditions, substrate accumulation was linear 
for at least 3 h (Fig. 1C) and then reached a plateau after 
overnight incubation (not shown). 

We next characterized several parameters of the FCoV 
nspl6-mediated MTase activity. An analysis of the pH depen- 
dence of the enzyme produced a bell-shaped profile, with max- 
imum activity observed at pH 7.5 (Fig. 2A). We also tested the 
dependence of the MTase on Mg?* and Mn?™, since these 
divalent cations previously were shown to promote O-methyl- 
ation by some MTases (42). FCoV nsp16 MTase activity was 
modulated by the presence of either of the cations Mn** (not 
shown) or Mg?* (Fig. 2B), with an optimum around a concen- 
tration of 1 mM. Finally, we tested the influence of the chain 
length of the capped substrate on nsp16 MTase activity. The 
protein was incubated with a variety of substrates, which dif- 
fered in the number of cytidines at their 3’ ends ("“°Gppp 
AC,,, where n = 1 to 7), and [*H]AdoMet for 4 h at 30°C. The 
nspl6-mediated methylation of RNA substrates again was 
measured using a filter binding assay. Figure 2C reveals a 
correlation between substrate length and MTase activity. The 
MTase activity increased with substrate length between 2 and 
5 nucleotides (n = 1 to 4) and then reached a plateau. In order 
to exclude that the observed differences are due to the incom- 
plete retention of small substrates on the filters used, the data 
were verified using HPLC analysis. To this end, nsp16 was 
incubated overnight with RNA and AdoMet, and the reaction 
products were separated by reverse-phase HPLC as described 
previously (58). The results were similar to those shown in Fig. 
2C (data not shown). We conclude that, under the experimen- 
tal conditions employed here, FCoV nsp16 exhibits the highest 
MTase activity at pH 7.5, in the presence of 1 mM Mg?*, on 
capped RNA substrates of the ”“°GpppAC,, type, where n 
equals at least four. 

FCoV nsp16 targets the 2'O position of the first transcribed 
nucleotide to produce a cap-1 structure. In order to identify 
the position(s) at which FCoV nsp16 methylates the 7“°Gpp 
pAC,; RNA substrate, we analyzed the reaction products by 
reverse-phase HPLC as described previously (58). As a con- 
trol, we included NSSMTasepy, which has been shown by mass 
spectrometry to methylate this substrate at the 2’O position of 
the first nucleotide (adenosine) and not onto downstream nu- 
cleotides (58). The retention time of the 7“°GpppAC, sub- 
strate was 29.7 to 29.8 min (Fig. 3A). Two additional peaks 
appeared when the reaction products of FCoV nsp16 were 
analyzed: one was clearly ahead of the AdoHcy reaction prod- 
uct, probably being a degradation product of the latter, and 
the other had a retention time of 33.9 min (Fig. 3C). These 
two peaks, but not the one corresponding to the substrate, 
also were observed among the products generated by 
NSSMTasepy, indicating that this enzyme methylated its sub- 
strate completely to yield ™“°GpppAzomeCs (58), which 
eluted at 33.9 min (Fig. 3B). These results indicated that both 
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FIG. 2. Biochemical characterization of FCoV nsp16 AdoMet-de- 
pendent MTase activity. The FCoV nsp16 MTase activity was mea- 
sured during a 120-min time course experiment by counting the 
amount of [*H]methyl transferred onto the RNA substrate (7“°Gpp 
pAC,). (A) Enzyme activity in reactions in Tris buffer (pH 7 to 9; 
triangles) or Bis-Tris buffer (pH 5 to 8; circles). (B) FCoV nsp16 
MTase activity in Tris-HCl buffer (pH 7.5) in the presence of increas- 
ing MgCl, concentrations. (C) FCoV nsp16 MTase activity after a 4-h 
incubation period using RNA substrates of increasing length (™“°Gp 
ppAC,.7). The cpm were normalized (100% corresponds to ™“°Gpp 
pAC,), and standard deviations were calculated from four indepen- 
dent experiments. 


enzymes catalyzed strictly the same reaction on this capped 
RNA, albeit with different efficiencies. 

The study of the cap structure was further refined as follows. 
We collected the 33.9-min peak material, digested it enzymat- 
ically to release the corresponding nucleosides, and analyzed 
the resulting mixture using HPLC (Fig. 3D, lower chromato- 
gram). The comparison to standard nucleosides (upper chro- 
matogram) indicated that all adenosines present in the reac- 
tion product were methylated at the 2'O position. Two 
additional peaks (3.7 and 7.4 min) also were observed, with the 
major one (3.7 min) corresponding to unmethylated C and the 
minor one (7.4 min) remaining ambiguous, corresponding to 


ysanB Aq Glog ‘G} yoreyy uO /Bio‘wse'IAl//:d}}Yy Wod papeojumoqg 


8076 DECROLY ET AL. 


™eGpppAC, 


29.8 min 


0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 
Minutes 


B  ™&GpppAC;+ NSSMTasepy 

™SGPPPAyoy4eC9 
Proteic material 33.9 min 
+ 


Remaining AdoMet 


0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 
Minutes 


C ™GpppAC; + FCoV nsp16 
0.32 


J. VIROL. 


™GPPPAyomeCs 
33.9 min 
™eG ppp AC, 
Proteic material 29.7 min 
ae 


Remaining AdoMet 


10 12 14 16 18 20 22 24 26 28 30 32 B 
Minutes 


T™MeG 
Cc ct 


2’OMe 


3.7 min Ani 


0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 
Minutes 


FIG. 3. FCoV nsp16 and NSSMTasepy methylate the 2'O position of the first transcribed nucleotide of the RNA. FCoV nsp16 or 
NSSMTasepy was incubated with ”“°GpppAC; overnight. The crude reaction mixture was analyzed using reverse-phase HPLC. The first section 
(in gray) indicates the removal of proteic material and remaining AdoMet by on-line cleaning on the precolumn (see Material and Methods). The 
gradient started after 5 min at 100% eluent A, with an increase to 10% eluent B after 25 min, reaching 30% after 45 min. (A) HPLC profile of 
’MeGpppAC;; (B) HPLC profile of 7“°GpppAC; incubated overnight with NSSMTasepy; and (C) HPLC profile of ”“°GpppAC, incubated 
overnight with FCoV nsp16. The peak eluting at 33.9 min was collected and digested with a mixture of nucleotide pyrophosphatase, phospho- 
diesterase I, and calf intestine phosphatase. (D) HPLC chromatogram of digestion products analyzed without the on-line cleaning procedure 
(lower chromatogram) compared to a mixture of standard compounds (upper chromatogram). The gradient started after 5 min at 100% eluent 
A, with an increase to 10% eluent B after 25 min and to 30% after 45 min. Absorbances (AU) were measured at 260 nm. 


7™eG or a mixture of ™°G and Cy,o,4.. We then performed a 
stoichiometric analysis of the cap components using the re- 
spective molar ratio corresponding to each nucleobase peak. 
The measured peak area (measured with the HPLC diode 
array detector) was divided by the molar extinction € at 260 nm 
that corresponded to each nucleoside (C = 7,200, G = 11,500, 
™eG = 10,100, and Aj.ome = 15,400). The relative molar 
ratios were found to be 0.82:1:6.2 for ™°G, Azone, and C, 
respectively. A cap-1 should yield ratios of 1:1:5, and a cap-2 in 
which ™°G comigrates with Cy5\. should yield 2:1:4. Taken 
together, these results demonstrate that FCoV nsp16 carries a 
2'O-MTase activity that is capable of converting a cap-0 to a 
cap-1 RNA structure. 

N7-methylated guanine of the cap structure is a binding 
determinant of FCoV nsp16. To gain insight into the mecha- 
nism of nspl6-mediated methylation, its binding to various 


substrates was analyzed. We conducted in vitro binding studies 
with small capped or uncapped RNAs of various lengths, with 
or without a methyl group at the N7 position. Using the T7 
DNA primase system (58), short a-**P-radiolabeled oligonu- 
cleotides were produced carrying ”“°Gppp, Gppp, or ppp at 
their 5’ ends. These RNAs were incubated with either FCoV 
nspl6 or NSSMTase,,, immobilized on Ni-NTA beads. After 
being incubated and washed, RNAs bound to the beads were 
resolved by PAGE and detected by autoradiography. As illus- 
trated in Fig. 4A, empty control beads bound little of the tested 
RNAs. In contrast, specific RNAs were retained by the differ- 
ent MTase-carrying Ni-NTA beads. NSSMTasep,, efficiently 
bound RNAs irrespective of their N7 methylation status 
(’“°GpppAC,., or GpppAC,.,) but did not bind uncapped 
pppAC,,, RNA. In contrast, FCoV nsp16 exclusively bound 
capped RNAs that had a methyl at the cap N7 position ("“°G 
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FIG. 4. Binding of noncapped RNAs or RNAs capped to FCoV 
nsp16 and NSSMTasepy. Short **P-radiolabeled RNAs were incu- 
bated with equal amounts of FCoV nsp16 MTase and NSSMTasepy, 
which previously had been immobilized on Sepharose beads. After 
three washes, the bound **P-labeled RNAs were separated by PAGE 
and detected by autoradiography. (A) PAGE analysis of pppAC,,, 
GpppAC,, and 7™°GpppAC,, bound to FCoV nsp16 MTase, 
NSSMTasepy, and empty Ni beads, which were used as negative con- 
trols (CTL—). The input (1/10) corresponded to 10% of the total 
radiolabeled RNA incubated with the MTase beads. (B) Binding of 
7™MeGpppAC,, to nsp16 MTase in the presence of increasing concen- 
trations of ™“°=GTP, ”“°GpppA, and Mg** cations. ND, not done, as 
the 20 mM point cannot be easily reached with commercially available 
™MEGpppA. 
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pppAC;3.,). Binding depended on RNA size, since capped 
RNAs with one or two 3’-terminal cytidines ("“°GpppAC, 5) 
did not bind efficiently to FCoV nsp16. Thus, requirements for 
efficient binding to nsp16 matched those favoring the MTase 
activity of the enzyme (compare Fig. 1B and 2C). We conclude 
that FCoV nsp16 specifically binds cap-0-bearing RNAs of at 
least 4 nucleotides in length. Both the RNA chain and the 
N7-methyl cap contribute to binding, but the presence of the 
latter is a prerequisite for binding. This substrate specificity is 
reminiscent of that of the VV 2'O-MTase VP39 (48) but quite 
different from that of the DV NSSMTase, which recognizes 
methylated and unmethylated capped RNA equally well (Fig. 
4A and reference 19) and harbors both N7- and 2'O-MTase 
activities in a single peptide chain (66). 

Since these results pinpointed nsp16 as a specific cap-0 bind- 
ing protein, we sought to determine whether ”“°GTP, GTP, or 
7MeGpppA could interfere with the binding of capped RNA, as 
previously shown for NSSMTasepy (20). To this end, we in- 
cubated nsp16-Ni-NTA beads with ”“°GpppAC;.¢ and mea- 
sured their binding capacity in the presence of an increasing 
concentration of GTP, ”™“°GTP, or ™*°GpppA. Figure 4B 
shows that neither GTP nor 7“°GTP (tested at up to 20 mM) 
had a clear effect on the interaction between nsp16 and RNA. 
The binding was slightly inhibited by 7“°GpppA, starting at 10 
mM. The binding of ™“°GpppAC;_¢ also was not affected by an 
increasing concentration of Mg?*, which was shown to stimu- 
late nsp16 MTase activity (Fig. 2B). Thus, Mg** ions were not 
essential for substrate recognition, and ““°GTP and ™°GpppA 
were poor competitors for binding. The latter property was 
consistent with the fact that ”“°GpppA could not efficiently 
serve as a substrate for nspl6. The VV VP39 2'O-MTase 
activity was not inhibited by 7“°GTP or ™“°GpppA (6). Inter- 
estingly, our nsp16 data contrast with the previously reported 
inhibition of the binding of NSSMTasepy to 7“°GpppAC;_, or 
GpppAC;_, by GTP and ™“°GTP (20). We conclude that the 
substrate binding site of the FCoV nsp16 MTase has a complex 
organization and interacts specifically with both cap-0 and nu- 
cleotides downstream of the RNA cap structure. 

Inhibition of FCoV nsp16 2'O-MTase activity by AdoMet 
and GTP analogues. Since viral capping enzymes are interest- 
ing targets for antiviral therapy (61, 85, 86), we screened po- 
tential inhibitors of FCoV nsp16 MTase activity, starting with 
AdoMet/AdoHcy and GTP analogues, which previously were 
identified as inhibitors of AdoMet-dependent MTases (45, 51, 
59-61). We first addressed whether GTP or 7™“°GTP could 
inhibit 2'O-MTase activity. For this purpose, nsp16 was incu- 
bated with 7“°GpppAC; and [*H]AdoMet in the presence of 
increasing 7“°GTP or GTP concentrations, and the incorpo- 
ration of label in methylated RNA was measured using the 
filter binding assay. Surprisingly, although GTP and 7“°GTP 
were unable to block the binding of nsp16 to its substrate even 
at 20 mM (Fig. 4B), we observed that both nucleotides inhib- 
ited MTase activity significantly, yielding IC, values of 1.51 + 
0.11 and 1.25 + 0.12 mM, respectively (Fig. 5A). This obser- 
vation suggests that these nucleotides do not act purely as 
competitive inhibitors of substrate binding. They might inter- 
fere with the binding of a specific part of the RNA substrate 
while leaving the overall binding unchanged. 

Using the same assay, we next tested whether FCoV nsp16 
was inhibited by a fixed (100 1M) concentration of different 
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GTP analogues, AdoMet, and the AdoHcy by-product of the 
methylation reaction (Fig. 5B). All tested compounds, includ- 
ing AdoHcy, which is known to be a potent inhibitor of the 
MTase family (6, 61), proved to be poor inhibitors of the 
MTase activity. The IC5, of AdoHcy was 144 uM (Fig. 5C), 
approximately 100-fold higher than the concentrations com- 
monly required to inhibit other MTases (6, 61). Likewise, sine- 
fungin, which previously was shown to be a potent inhibitor of 
the MTases of Newcastle disease virus and VV (IC; values of 
150 and 75 nM, respectively [60]), was found to be a very poor 
inhibitor for FCoV nsp16. Thus, we conclude that nsp16, un- 
like its distant viral homologues, can efficiently discriminate 
AdoMet and its analogues, suggesting that although various 
AdoMet-dependent MTases share a similar structural organi- 
zation (see below), nsp16 evolved a specific variant of the 
conserved AdoMet binding site. 

Identification of key residues for FCoV nsp16 2’O-MTase 
activity and substrate binding. To identify residues involved in 
catalysis and substrate binding, we engineered and character- 
ized a set of FCoV nsp16 point mutants. The positions to be 
mutated were selected using a comparative analysis of nsp16 
and other, better characterized MTase homologues. We 
started from a structure-based sequence alignment of four 
RNA 2'O-MTases with known tertiary structures. These en- 
zymes of viral and cellular origin were aligned with nsp16 of 
SARS-CoV and FCoV, and for the latter a secondary structure 
prediction was generated (Fig. 6A). Due to the considerable 
divergence of the protein sequences used, only limited conser- 
vation is evident in the alignment, in particular in regions 
corresponding to the MTase core, which is structurally con- 
served among AdoMet-dependent MTases (14, 52). In line 
with the original analysis of CoV nsp16 (76, 82), four out of six 
invariant residues in the alignment correspond to the K-D-K-E 
catalytic tetrad of RNA 2'O-MTases (12, 19). These four res- 
idues (K45, D129, K169, and E202 in FCoV) were selected for 
mutagenesis. To identify other potentially important residues, 
the tertiary organization of the FCoV nsp16 MTase core was 
modeled through different Web servers (see Materials and 
Methods). Although the VP39 and FCoV nsp16 MTases share 
the same binding specificity for cap-0 structures, they were 
found to be too divergent for any server to produce a mean- 
ingful nsp16 model. The closest template proposed by several 
servers was the rRNA 2'O-MTase FtsJ (11), which therefore 
was used to model FCoV nsp16 using HHpred (77). An in- 
spection of the final model (Fig. 6B) of the MTase core of 
nspl6 (residues K45 to G221) and other MTase structures 
(Fig. 6C and D) showed that in all structures but that of VP39, 
the spatial equivalent of FCoV nsp16 residue D113 (conserved 
in CoV nsp16) is also an Asp residue (11, 19), which is located 
close to the AdoMet. Therefore, residue D113 of FCoV nsp16 
was selected for mutagenesis. 

The identification of residues that might be involved in sub- 
strate binding was less straightforward. Our starting point was 
the observation that VV VP39 shares with nsp16 a high spec- 
ificity for the cap-0 structure. In the structure of VP39 com- 
plexed with a 7™“°G-capped RNA (Fig. 6C) (33), the methyl- 
ated base is stacked between two aromatic side chains (Y22 
and F180) (Fig. 6A and C). Moreover, in VP39 the methyl 
group is in contact with residue Y204 (Fig. 6A and C) via van 
der Waals interactions (33, 64). Finally, the carboxyl groups of 
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FIG. 5. Inhibition of FCoV nsp16 MTase activity by GIP and 
AdoMet analogues. The nsp16 MTase was incubated with [*H]AdoMet 
and ™°GpppAC, in the presence of various potential inhibitors. MTase 
activity was determined by measuring the radioactivity associated with the 
RNA substrates using filter binding assays. (A) Inhibition curve of nsp16 
activity by GTP (O) and ™“°GTP (A). The nsp16 activity obtained in the 
absence of inhibitors was set to 100%. The standard deviations from three 
independent experiments are shown, and the IC;, value was calculated as 
indicated in Materials and Methods. (B) Inhibition by various AdoMet 
analogues at a concentration of 100 4M. Activity of 100% corresponds to 
the nsp16 activity on the 7”“°GpppAC, substrate in the absence of inhi- 
bition. CTL, control. SIBA, 5’-deoxy-5’-S-isobutylthioadenosine; TP, 
triphosphate. (C) Inhibition curve of FCoV nsp16 activity by the reaction 
product (AdoHcy). The ICs, values were determined as described for 
panel A. 
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residues D182 and E233 form hydrogen bonds to the NH and 
NH2 of the guanosine in VP39. Most of these residues are 
located in regions for which the counterparts in the FCoV 
nsp1l6 could not be modeled by HHpred. The only exception 
was the loop region immediately downstream of B-strand 6, in 
which we selected residues F173 and W175 for mutagenesis, as 
they are a putative counterpart of VP39 F180 (Fig. 6B). To 
probe for cap binding residues in the most divergent N- and 
C-terminal regions, we limited our selection to residues con- 
served in CoV that are physicochemically similar to those in- 
volved in substrate recognition in VP39. Aromatic residues 
W4, Y14, and Y29 in the nsp16 N terminus were selected as 
putative equivalents of Y22 and Y204 in VP39 and residues 
D221 and D247 in the nsp16 C terminus (Fig. 6A) as candi- 
dates for the role of E233 in VP39. 

Of the 12 nsp16 mutants engineered, all but W175A were 
expressed at a level sufficient to allow the purification of the 
enzyme by our two-step procedure. Consequently, W175A 
could not be analyzed further. Both their 7“°GpppAC;, bind- 
ing capacity and their MTase activity on capped RNA oligo- 
nucleotides ("“°GpppAC; and GpppAC;) were compared to 
those of the wt control (Fig. 7). For none of the mutants was 
there a change in the dependence on a N7-methylated cap or 
the selectivity compared to that of an unmethylated cap (Fig. 
7A). The MTase activity was abolished or severely reduced by 
all mutations affecting the putative K-D-K-E catalytic tetrad, 
supporting their structure-based assignment as being function- 
ally essential (Fig. 6A). One of these mutations (K45A) also 
abolished the binding of capped RNA (Fig. 7C), while muta- 
tions at three other positions reduced it to 20 to 50% of the 
level observed for the wt protein. This result indicated that 
K45, and possibly other catalytic residues, contribute to the 
specific recognition of cap-0 RNA. In the case of D113, which 
may be part of the predicted AdoMet binding pocket, replace- 
ment by Ala reduced 2’O-MTase activity to about 10% of the 
wt level without an apparent effect on substrate binding. In 
contrast, Ala replacements of the aromatic residues located in 
the N-terminal subdomain (Y14 and Y29) and downstream of 
B-strand 6 (F173) also reduced the binding of nsp16 to its 
substrate to 10 to 20% of the wt level (Fig. 7C). The effects on 
the MTase activity varied considerably, from complete inhibi- 
tion (Y29A) to a modest 50% reduction (F173A) or a minor 
~18% drop (Y14A). Mutations W4A, D221A, and D247A 
formed a group in which the effect of replacement by Ala on 
both RNA binding and 2'’O-MTase activity was comparable, 
ranging from no change to 55% inhibition at most. 

We conclude that the FCoV nsp16 2’O-MTase employs the 
canonical K-D-K-E catalytic tetrad that consists of K45, D129, 
K169, and E202. Residue D113 also may contribute to catalysis 
but not to RNA binding. Two residues (K45 and Y29) for 
which mutagenesis severely compromised both MTase activity 
and substrate binding may be involved in both substrate rec- 
ognition and catalysis or, alternatively, only in the specific 
recognition of cap-0 RNA. 


DISCUSSION 


Upon CoV infection, replicases ORFla and ORF 1b of the 
positive-strand RNA genome are translated, presumably by 
using a cap-dependent mechanism (44, 80). Although the cap 
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is expected to be essential for ribosome binding and efficient 
viral mRNA translation (24, 75), the RNA capping mechanism 
of CoV is essentially uncharacterized. It supposedly involves 
four enzymatic activities (24), including an RTPase (40, 41), a 
GTase, and two MTase activities. 

The ORF1b-encoded replicase subunits of CoV (nsp12 to 
nspl6) include five RNA-processing enzymes that have 
uniquely segregated in the genomes of CoV and other large 
nidoviruses (28, 76). In recent years, the activities of four of 
these enzymes, but not the nsp16 MTase, were characterized 
for one or more CoV. In this report, we document the MTase 
activity of nsp16, the C-terminal subunit of the CoV replicase. 
FCoV nsp16 specifically and selectively binds capped RNAs of 
3 to 6 nucleotides in length that carry a methyl group at the N7 
position of the guanosine cap, a chemical entity known as the 
cap-0 structure (’“°GpppAC;_,). FCoV nsp16 catalyzes the 
transfer of a methyl group from the AdoMet donor to short 
capped RNAs, but only when they carry a methyl group at the 
N7-guanine position. 

Different 2'O-MTases were reported to target either the first 
(71, 72) or subsequent transcribed nucleotides of capped 
RNAs (32). We have identified the methylation position using 
an enzymatic digestion of the reaction products in conjunction 
with HPLC analysis. The latter demonstrated that the nsp16 
MTase targets the 2’O position of the first nucleotide, thus 
converting ™“°GpppAC,, into ”“°GpppA3-omeC,, Thus, FCoV 
nsp16 belongs to the family of AdoMet-dependent mRNA cap 
2'O-MTases that specifically methylates the first endogenous 
nucleotide of the cap-0 structure to produce a cap-1 structure. 
These findings support the proposed involvement of nsp16 in 
CoV translational control (76, 82) and the suggestion that CoV 
mRNAs carry a cap structure (44, 80). A more extensive anal- 
ysis of nsp16 is required to verify whether this enzyme also 
methylates other nucleotides in virus and/or cellular RNAs to 
regulate other processes in cooperation with virus RNA-pro- 
cessing enzymes, as proposed previously (76). With these find- 
ings, two of the four cap-forming enzyme functions now have 
been reported for CoV: an RTPase activity was mapped to the 
multifunctional nsp13 helicase protein of HCoV 229E (41) and 
SARS-CoV (40) and the 2’O-MTase activity reported here for 
FCoV nsp16. 

Our finding that the FCoV nsp16 2’O-MTase strongly pre- 
fers substrates carrying a cap-0 structure over those that have 
a nonmethylated cap indicates that cap maturation in CoV 
follows the canonical order of methylation steps that have been 
described for other biological systems. This will have to be 
compared to the cap methylation recently documented for 
flaviviruses, in which the NSS MTase recognizes both methyl- 
ated and nonmethylated cap structures (this work and refer- 
ence 20) and catalyzes the methylation of both the N7 position 
of the guanine as well as the 2’O position of the nucleotide 
ribose (66, 86). It is possible that nsp16 also possesses N7- 
guanine MTase activity in the presence of longer substrates, 
e.g., RNA sequences and/or structures specific for the 5’-un- 
translated region of the viral genome. For example, the spec- 
ificity of the West Nile virus (WNV) NSSMTase for WNV- 
specific RNA sequences recently was reported (16, 66). If such 
a dependence exists in the case of CoV, RNA sequences down- 
stream of the cap structure might represent regulatory ele- 
ments inducing specific N7 or 2'O-methylation events. 
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FIG. 6. Comparative analysis of viral MTases. (A) Structure-based alignment of the MTase core domain of FCoV nsp16 with SARS-CoV nsp16 
(82), the VV mRNA MTase VP39 (34), NSSMTasepy (19), the rRNA MTase FtsJ (11), the C-terminal rRNA MTase domain of Mj0697 (83), and 
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Residues potentially involved in substrate binding and the 
methyl transfer reaction were probed by mutagenesis. In ad- 
dition to the catalytic tetrad K-D-K-E, residue D113 of the 
AdoMet binding site and residue Y29 are essential for 2’O- 
MTase activity. Similar results were reported recently for the 
WNV MtTase. Using alanine-scanning mutagenesis, it was 
shown that its K-D,,,-K-E catalytic tetrad is essential for 2'O 
methylation activity, whereas the N7-MTase activity depends 
mainly on D,,, (86). The catalytic residues are adequately 
positioned in the three-dimensional structural model of FCoV 
nsp16 (Fig. 6B). Interestingly, we noted that some residues of 
the catalytic tetrad also are important for the efficient binding 
of the capped RNA substrate. This observation suggests that 
the binding of cap RNA involves multiple contact regions. 
Accordingly, binding assays revealed that a shorter cap-1-con- 
taining RNA, ”™“°GpppAC,.», did not bind to FCoV nsp16, 
whereas longer molecules did. The more efficient binding of 
longer RNAs suggests that the presence of nucleotides down- 
stream from the cap binding site is necessary to stabilize the 
interaction between FCoV nsp16 and its RNA substrate. This 
also was confirmed by the poor inhibition of binding observed 
in the presence of increasing concentrations of 7“°GpppA cap 
analogue and ’“°GTP. In particular, the K45A substitution 
was quite spectacular in terms of its effect on both MTase 
activity and RNA binding. A single nsp16 point mutation abol- 
ishing RNA binding is remarkable. Indeed, ”™“°GpppAC; 
should interact with many residues of the enzyme, each con- 
tributing moderately to RNA binding. Could K45 be an essen- 
tial residue responsible for specific cap recognition? The su- 
perimposition of the VP39 and nsp16 models shows that the 
methylated guanosine cap binding site is remote from K45 
(Fig. 6B and C), suggesting that K45A has a drastic indirect 
effect on the putative RNA binding groove but not a direct 
effect on the binding of the methylated cap. In contrast, the 
fact that D113A did not influence substrate binding is consis- 
tent with its specific position within the putative AdoMet- 
binding site. The corresponding residue in NSSMTasepy is 
located on the side of the AdoMet cofactor, not adjacent to the 
positively charged, putative RNA-binding groove (20). 

Our mutagenesis study confirmed a contribution of at least 
one aromatic amino acid to the specificity of methylated cap 
recognition. Since Y29A abolished both cap-0 recognition and 
the 2'O-MTase activity of FCoV nsp16, we propose that this 
residue plays an essential role in specific cap recognition 
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through a stacking interaction with cap-O RNA, which in turn 
might allow the correct positioning of capped RNA into the 
catalytic site. Unfortunately, Y29 could not be visualized in the 
structural model. In contrast, the Y14A and F173A replace- 
ments did not reduce 2’O-MTase activity as extensively as they 
inhibited RNA binding. This observation may be explained by 
the stringent washing conditions in our binding assay and sug- 
gests that both residues contribute to cap recognition. The cap 
binding site of the VV VP39 MTase previously was shown to 
depend on two aromatic side chains (Y22 and F180) that are 
necessary to recognize the cap-0 structure by stacking interac- 
tions (35), and the replacement of F180 drastically inhibited 
2'O-MTase activity (36). If we assume that a cap binding site is 
located in the same area in both FCoV nsp16 and VV VP39, 
F173 is adequately positioned in the three-dimensional model 
of nsp16 (Fig. 6B) to bind to the cap-0 structure. Therefore, we 
hypothesize that FCoV nsp16 stacks the methylated guanine 
moiety of the cap-0 structure between two aromatic residues, 
similarly to what was found for VV VP39 and the mammalian 
cap binding translation initiation factor eIF4E (32, 35, 64). As 
the most likely candidates for this role, we propose residues 
Y29 and F173 without formally excluding Y14. 
Virus-encoded RNA cap 2'O-MTase activities previously 
have been identified for various other virus genera, such as 
poxviruses (VP39 [71, 72]), the double-stranded reoviruses (50, 
65), minus-strand ssRNA viruses such as VSV (47), and posi- 
tive-strand ssRNA viruses like flaviviruses (NSSMTase [19, 58, 
86]). The exact role of mRNA cap 2’O-methylation in the virus 
life cycle is still an open question. It is likely that this function 
confers a replicative advantage on the virus, since 2’O-MTases 
are widely conserved among viruses that have cap-bearing ge- 
nomes and/or mRNAs. For the WNV MTase, it has been 
shown that mutations abolishing the N7 and/or 2'O-MTase 
activity have a detrimental effect on the replication of a lucif- 
erase-expressing RNA replicon (66). In contrast, when viral 
replication was tested directly, a substitution inactivating the 
N7-MTase activity substantially attenuated replication, 
whereas 2'O-MTase knockouts had more moderate effects (66, 
86). These observations suggest that compounds specifically 
inhibiting cap-methylating enzymes, either N7-MTase or 2'O- 
MTase or both, could act as potent antiviral agents. Accord- 
ingly, a series of AdoMet analogues containing modifications 
of the nucleoside moiety showed potent and selective antiviral 
activity against herpes simplex virus, VV, and VSV (5, 61). In 


the mRNA MTase domain I of reovirus (Reov) \2 (67). The alignment was generated manually using Seaview (25) and was presented with the 
ESPript program (30). The positions of conserved residues are highlighted by a red background, and the regions outlined in blue indicate high 
sequence similarity (>70%). The secondary-structure elements of the VP39 and NSSMtasepy structures and their names are given below the 
alignment. The predicted FCoV nsp16 secondary structure generated by the PSIPRED protein structure prediction server (10) is indicated above 
the alignment. The position of the mutations characterized in this study are indicated by triangles above the FCoV nsp16 sequence: salmon 
triangles for the putative cap-0 binding site, light blue/turquoise for the catalytic site, and yellow for the AdoMet binding site. The conserved 
catalytic active-site residues K-D-K-E of RNA 2'O-MTases are annotated by blue stars, and the VP39 residues participating in cap recognition are 
indicated with salmon circles. Also shown are the structural comparative analyses of viral Mtases. (B) Ribbon representation of the FCoV nsp16 
model as generated on the basis of the structure of rRNA 2’O-MTase FtsJ (see Materials and Methods). AdoMet was placed by the superim- 
position of the FCoV nsp16 model and the structure of rRNA 2’O-MTase FtsJ in complex with the cofactor. Structures of VP39 (PDB ID 1AV6) 
(C) and NSSMTasepy (PDB ID 1L9K) (D) also are shown. A yellow stick representation is used to show AdoMet and the reaction product 
AdoHcy, with oxygen in red and nitrogen in blue. The side chains and carbon backbone of conserved catalytic residues K-D-K-E are represented 
by sticks in blue/turquoise (K45 was manually added in the FCoV nsp16 model). The lateral chain of nsp16 D113, which is involved in AdoMet 
binding, is shown in yellow. The VP39 side chains of residues involved in the binding of the ™“°G cap are represented in salmon, and the cap 
analogues are in light gray. 
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FIG. 7. Alanine-scanning mutagenesis of the FCoV nsp16. Activity 
and substrate binding of FCoV nsp16 mutants. (A) Equal amounts of 
the different nsp16 mutants were incubated with GpppAC; and ™“°G 
pppAC; in the presence of [*H]AdoMet. The methyl transfer to the 
RNA substrate was monitored during a 240-min time course experi- 
ment and detected using a filter binding assay that measured the 
amount of radioactivity, in cpm, that was transferred to the RNA 
substrates. MTase activity of 100% was arbitrarily attributed to the wt 
nsp16 activity when using the “°GpppAC; substrate. The bar graph 
presents the results of three independent experiments. (B) The puri- 
fied nsp16 mutants after separation on SDS-12% Page and Coomassie 
blue staining. (C) The binding of ’“°GpppAC; to nsp16 mutants was 
determined as described in the legend to Fig. 4. Binding data were 
normalized (100% corresponds to ’“°GpppAC; bound to wt nsp16). 


this study, we tested whether some of these AdoMet analogues 
(e.g., sinefungin, which has an IC;, of 39 nM for VV VP39 
[60]) would block nsp16 MTase activity, but all of the mole- 
cules tested proved to be weak inhibitors. 

Our observations suggest that the CoV nspl6 MTase is a 
unique enzyme: it shares its structural organization and spec- 
ificity for the cap-0 structure with VV VP39, but the divergence 
of the structurally conserved AdoHcy binding site results in 
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unusually weak inhibition profiles for AdoMet analogues. Fur- 
thermore, nsp16 has evolved a unique substrate binding site 
that cooperates with the canonical catalytic tetrad that is con- 
served in other 2’O-MTases to recruit its substrate. Finally, 
despite the inability of GTP and ™“°GTP to block cap-0 bind- 
ing, both nucleotides inhibit FCoV nsp16 activity at similar 
concentrations, suggesting that the nucleotide binding site is 
distinct from the site binding the methylated cap. Further func- 
tional and structural exploration of FCoV nsp16 will be nec- 
essary to define this putative binding site, which is of potential 
interest for inhibitor design. 
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